Dynamic Entropy as a Measure of Caging and Persistent Particle Motion in 
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The length-scale dependence of the dynamic entropy is studied in a molecular dynamics simulation 
of a binary Lennard- Jones liquid above the mode-coupling critical temperature T c . A number of 
methods exist for estimating the entropy of dynamical systems and we utilize an approximation 
based on calculating the mean first-passage time (MFPT) for particle displacement because of its 
tractability and its accessibility in real and simulation measurements. The MFPT dynamic entropy 
S(e) is defined to equal the inverse of the average first-passage time for a particle to exit a sphere of 
radius e. This measure of the degree of chaotic motion allows us to identify characteristic time and 
space scales and to quantify the increasingly correlated particle motion and intermittency occurring 
in supercooled liquids. In particular, we identify a "cage" size defining the scale at which the 
particles are transiently localized, and we observe persistent particle motion at intermediate length 
scales beyond the scale where caging occurs. Furthermore, we find that the dynamic entropy at 
the scale of one interparticle spacing extrapolates to zero as the mode-coupling temperature T c is 
approached. 



T3 

C 

O 
o 



> 

m 
m 
o 



i 

-a 
c 

o 
o 



X 



I. INTRODUCTION 

It has been suggested that the glass transition in cooled 
liquids is a dynamic transition from an ergodic to a non- 
ergodic state. For example, the ideal mode-coupling the- 
ory predicts that the molecules of simple liquids become 
increasingly "caged" by surrounding molecules, result- 
ing in an ergodic to non-ergodic transition at a criti- 
cal temperature T c at which the fluid molecules become 
permanently localized (i.e., the self-diffusion coefficient 
vanishes) jjj. Although a tendency toward particle lo- 
calization for increasingly long times has been observed 
in simulations and experiments on supercooled liquids, 
particle localization and structural arrest does not actu- 
ally occur at the extrapolated temperature T c because 
the particles are eventually able to "escape" their cages. 
Recent simulations have also shown the tendency for par- 
ticle motion to occur in an increasingly correlated way in 
supercooled liquids a feature emphasized by the 

older, phenomenological Adam-Gibbs model of glass for- 
mation P]. The observed greater particle mobility near 
T c is presumably a consequence of the increased collec- 
tive motion of cooled liquids ( "hopping" in the extended 
version of the mode-coupling theory jl|) which restores 
the ergodicity of the liquid for some temperature range 
below T c . This thermally activated collective motion ap- 
parently postpones the ergodic to nonergodic transition 
to a lower temperature. In the Adam-Gibbs model 
this lower temperature corresponds to the conjectured 
"ideal" glass transition temperature T Q , where the equi- 
librium configurational entropy extrapolates to zero |pL0|| . 

If glass formation indeed represents an ergodic to non- 
ergodic dynamic transition, then it is important to de- 



fine a dynamical measure of order that quantifies both 
the "closeness" of the transition pi]] , and the degree of 
correlated motion in an equilibrium glass-forming liquid. 
Ergodic theory provides us with a natural measure in the 
form of the "dynamic entropy" p^TTt . 

The concept of dynamic entropy was introduced by 
Shannon in his theory describing the capacity of ideal 
communication devices to transmit information [p"8| . 
This idea was later developed by Kolmogorov and others 
|]l3| into a general measure of the "degree of random- 
ness" or "degree of chaos" of dynamical systems. Ac- 
cording to Pessin's theorem p9[ , the Kolmogorov-Sinai 
dynamic entropy \iks for & Hamiltonian dynamical sys- 
tem equals the sum of the positive Lyapunov exponents 
p9|^0[ . These exponents are measures of the "instabil- 
ity" of the system evolution p^pl. Dynamic entropy 
extends the equilibrium definition of entropy from sta- 
tistical mechanics to the time domain. The dynamical 
entropy provides an estimate of the rate of growth of 
"information" (per unit time) required to describe the 
evolution of a dynamical system |l4||l5|,^l| and it is also 
a measure of the "complexity" of a dynamical system 
p2| . The dynamic entropy characteristically decreases 
as a system orders and its exploration of its phase space 
becomes more restricted |2^,Q. Thus, the dynamic en- 
tropy decreases as a fluid crystallizes or a spin system 
orders |§4|-§7||. 

The Kolmogorov-Sinai dynamic entropy has some 
shortcomings in the description of complex configura- 
tional changes that occur in supercooled liquids. In par- 
ticular, hxs diverges for the ideal process of Brownian 
motion (due to the non-differentiability of the trajecto- 
ries) p| , p8p9| ]. Consequently, we must anticipate diffi- 
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culties in applying dynamic entropy to quantify particle 
motions at large length and time scales in the case of 
supercooled liquids. Recently, there has been an impor- 
tant generalization of the dynamic entropy concept that 
provides a "bridge" between microscopic dynamical sys- 
tem descriptions and macroscopic stochastic descriptions 
of liquid dynamics. This generalization recognizes that 
the amount of information required to describe the paths 
of a stochastic process depends strongly on the length 
scale of observation e. The e-dependent dynamic entropy 
h(e) of Gaspard and Wang and others (l4| (also called 
"e-entropy") reduces to the hxs entropy in the limit of 
small e, 



lim h(e) 

£-+0 



(1.1) 



and is well-defined for idealized stochastic processes at a 
fixed, nonvanishing e. The dynamic entropy of a Brown- 
ian particle /is(e) obeys the scaling relation 



h B (e) oc e' 



e > 0, 



(1.2) 



where the proportionality constant is fixed by the par- 
ticle diffusion coefficient jL4|. As mentioned above, h(e) 
for stochastic particle motion diverges as e — > and the 
exponent reflects the fractal dimension of the particl e tr a- 
jectories |l^ , |30|l . Specifically, the exponent 2 in Eq. F2 is 
the Hausdorff dimension of a Brownian path in three di- 
mensions |}l| , and in the limit of perfectly coherent (bal- 
listic) particle motion this exponent is 1. In idealized 
stochastic processes (e.g. fractional Brownian motion, 
Levy flights, etc.) the exponent in Eq. 1.2 can be identi- 
fied with the path Hausdorff dimension [p^ , pr||30|j32| , and 
can take values intermediate between 1 and 2. This ex- 
ponent reflects the "degree of persistence" in the particle 
displacement relative to Brownian motion. 

The scale dependent dynamic entropy h(e) for com- 
plex dynamical systems such as liquids depends strongly 
on the observational scale e. At very small e the mi- 
croscopic chaotic motion of the molecules is observed, 
so that h(e) varies slowly with e. The decorrelation of 
particle velocities in a liquid occurs at a time and space 
scale corresponding to the average interparticle "collision 
time" and h(e) starts varying with e as this decorrelation 
occurs. This helps us to identify a characteristic space 
and time scale over which the bare microscopic dynamics 
can be coarse-grained by a stochastic description. Cor- 
relations associated with particle displacement arise at 
longer times in cooled liquids and h(e) also helps us in 
determining the spatial and time scales over which these 
correlations occur. h(e) thus provides a measure of the 
degree of chaotic motion appropriate to the description 
of real systems at arbitrary observational scales and is an 
attractive tool for quantifying the increasingly restricted 
motion in cooled liquids. It is notable that its definition is 
not restricted to circumstances where statistical mechan- 
ical equilibrium exists, so that this measure of the degree 
of chaos extends to non-equilibrium situations such as 
the glass state and turbulent fluids Q . 



The calculation of h(e) O] (or Hks) is generally dif- 
ficult, especially in cases where h KS is small and long 
computational times are required for its accurate deter- 
mination @|||. In the present paper, we utilize a simple 
approximation for h(e) that has the advantage of being 
accessible in experiments on real materials and computer 
simulations |14B3[. Provided that the spatial scale e is 
not too small fll4|| , h(e) can be approximated by enclos- 
ing the particle position at time t = by a sphere of 
radius e centered on the particle, and then determining 
the time t at which the trajectory first arrives at the 
threshold distance e (see Fig. 1). We average this "first- 
passage time" over all particles in the liquid to obtain 
the mean first-passage time (MFPT) r(e) and we define 
the "MFPT dynamic entropy" 5(e) as, 



5(e) = l/r(e), 



where 



r(e) 



dtP e (t) t. 



(1.3) 



(1.4) 



P e (t)dt is the probability that the particle reaches the 
distance e between t and t + dt. The dynamic entropy 
5(e) is thus one measure of the average "escape rate" 
of a particle from its local environment [53. We note 
that although the definition of 5(e) is motivated by dy- 
namical systems theory concepts, this property defines 
an independently interesting measure of correlated mo- 
tion in liquids that does not rely on the approximation 
relating 5(e) to h(e). 




FIG. 1. Schematic of a particle trajectory in a cooled liq- 
uid. The solid line represents the e-sphere (colored gray). t(e) 
is the first-passage time for the particle to reach the sphere 
boundary. Filled circle denotes initial particle position, and 
open circle denotes particle at first-passage time. 

In this paper we utilize 5(e) to identify characteristic 
space and time scales in the particle dynamics and to 
quantify the increasingly correlated motion observed in 
previous analyses of the same simulations considered in 
the present paper These studies indicated the de- 

velopment of large scale dynamical heterogeneity and the 
nature of this dynamical heterogeneity has been exam- 
ined in a series of recent papers that are complementary 
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to the present work There it was established that 

transient clusters of highly "mobile" particles form in the 
cooled liquid and that the average size of these clusters 
grows rapidly as T c is approached pL A pair distribu- 
tion function for particle displacements was defined and 
this quantity exhibits a growing length scale upon cool- 
ing that reflects the clustering of mobile particles . 
The growing length scale is time-dependent and attains 
a peak value at a time in the a-relaxation regime ||,|7j . 

It has also been shown in the present liquid that the 
particles within the mobile particle clusters move in co- 
operatively rearranging "strings" Notably, the 
string- like collective motion also begins well above T c , 
but the strings themselves exhibit no tendency to grow 
rapidly near T c . Instead, the length distribution of the 
strings is found to be nearly exponential and a similarity 
of this distribution to that commonly observed in equi- 
librium polymerization has been noted Donati et al. 
have suggested that the growing barrier height to par- 
ticle motion is proportional to the average string length, 
which would imply that these string-like motions have a 
basic significance for understanding transport in cooled 
liquids. Thus, string-like correlated particle motion ap- 
pears to be an important mode of motion in our cooled 
liquid and part of the motivation of the present work 
is to better characterize the development of this type of 
collective motion. We are also interested in the extent 
to which particle displacement becomes intermittent in 
time in cooled liquids, since a growing intcrmittency in 
particle motion has been suggested to underlie the glass 
transition ]3^ , |37[ |. 

The paper is outlined as follows. In Section II we re- 
view some details of the MD simulation data utilized in 
this work. Section III examines 5(e) over a broad range 
of scales, and dynamical regimes are defined where the 
motion is ballistic, transiently localized, persistent and 
diffusive. These regimes are examined in separate sub- 
sections. We summarize our findings in Section IV. 

II. SIMULATION DETAILS 

The system studied is a three-dimensional binary mix- 
ture of 8000 Lennard-Jones (LJ) particles in which the 
sizes of the particles and the interaction parameters are 
chosen to prevent crystallization and demixing |3^] . The 
size of the A-particles is about 10% larger than that of 
the B-particles (while the mass is the same) and the par- 
ticles have a relative concentration 80:20 of A-particles 
to B-particles. We report our results in dimensionless 
LJ units p8|. The system was equilibrated at different 
temperatures T in the range (0.451,0.550). The density 
p varied from 1.09 particles per unit volume at the high- 
est temperature to 1.19 at the lowest T simulated. For 
reference, the mode-coupling temperature T c for this sys- 
tem is T c = 0.435 at p ~ 1.20 ||@J§, so all the simula- 
tion data analyzed here is well above the glass transition. 



Configurational histories for up to 4 ■ 10 6 molecular dy- 
namics time-steps following equilibration were stored for 
each run. Following equilibration in the NPT and NVT 
ensembles, the trajectories were calculated in an NVE 
ensemble, and snapshots containing the particle coordi- 
nates and velocities were taken at logarithmic time inter- 
vals during the run. In this stage, the equations of mo- 
tion were integrated using the velocity Verlet algorithm 
with a step size of 0.0015 at the highest temperature, and 
0.003 at all other temperatures. Adopting argon values 
for the LJ parameters of the large particles implies an 
observation time of « 26 ns for the coldest T. All data 
presented here are calculated for the majority (A) par- 
ticles only, except where otherwise noted jf0|. Further 
details of the simulation can be found in Rcf. pi . 




1Cf 2 1Cf 1 10° 10 1 10 2 10 3 10 4 
t 

FIG. 2. Mean square displacement of the majority species 
(A particles) vs. time for different T. 

Over the temperature-density regime studied, the sys- 
tem exhibits the usual features of a fragile EJ , glassform- 
ing liquid. For example, the mean square displacement 
<r 2 (i)} = (]^Eil kiW " r*(0)| 2 ) for the A particles 
is shown in Fig. || for different T. Here r^(t) is the po- 
sition of particle i at time t, Na is the number of A 
particles (6400), and (• • •) denotes an ensemble average. 
For each state point, a "plateau" exists in both the mean 
square displacement and the self-part of the intermedi- 
ate scattering function F s (q, t) as a function of t (see || ) . 
The plateau in Fig. || separates an early time "ballistic" 
regime from a late time diffusive regime. The plateau 
is interpreted as implying "caging" of the particles and 
this phenomenon is typical for liquids of low tempera- 
ture or high density. Over the range of T studied, the 
a-relaxation time r a , describing the decay of F s (q, t) (at 
the value of q corresponding to the first peak in the static 
structure factor), increases by 2.4 orders of magnitude, 
and follows a power law t q ~ (T — T c ) -7 , with T c ~ 0.435 
and 7 ~ 2.8. The simulated liquid states analyzed here 
therefore exhibit relaxation behavior characteristic of a 



3 



supercooled liquid. No long range structural correlations 
due to density or composition fluctuations are apparent 
in the simulation data [Si. 



III. CHARACTERISTIC DECORRELATION 
TIME AND SPACE SCALES 

In Fig. H we show the MFPT r(e) plot for six different 
runs corresponding to varying the temperature of the sys- 
tem from T = 0.550 to T = 0.45f . In the inset we show 
the dynamic entropy S(e) = l/r(e). Note that the vari- 
ation of r(e) with e exhibits similar qualitative trends to 
the variation of t with (r 2 (t)) shown in Fig. g. 

For small e, corresponding to the inertial regime, S(e) 
is insensitive to temperature. At intermediate e values 
we see a decrease of S(e) with decreasing T and an in- 
crease in the magnitude of the slope in the log- log plot. 
A strong temperature dependence of S(e) is apparent at 
a scale on the order of one interparticle distance (e = 1). 
On these larger scales, we show in a later subsection that 
S(e) exhibits a power-law scaling with e and reduced tem- 
perature. 
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FIG. 3. Mean first-passage time r of the majority species 
(A particles) versus e, for different T. Inset: Mean 
first-passage time dynamic entropy S(e). Compare the 5(e) 
variation observed here for a cooled liquid with the dynamic 
entropy h(e) calculated for a one-dimensional model map ex- 
hibiting diffusion at long times (see Fig. 25b of Gaspard and 
Wang [8]). 

It is apparent from Fig. || that the particle displace- 
ment in this cooled liquid is not Brownian over most 
of the simulation time scales, and it is conventional to 
quantify this deviation by a "non-Gaussian parameter" 
a(t) involving the moments (r 2 (t)) and (r (t)) of the 
self part of the van Hove correlation function G s (r, t) = 

(ft E^i S(r - (r;(i) - r 4 (0)))). The parameter a(t) is 
defined as, 



= 3(r 4 (t)) _ 
a{ > ~ 5(r 2 (t)) 2 ' 



(3.1) 



and vanishes for Brownian motion. a(t) is shown in Fig. [| 
together with r(e) for the coldest run (T = 0.451); note 
that for a(t), time t is plotted on the ordinate axis. This 
comparison allows us to identify four regimes: An inertial 
regime (Regime I) where the non-Gaussian parameter is 
small; a "localization" regime characterized by a large 
value for the slope in the r(e) log- log plot and by a grow- 
ing a(t) (Regime II); a regime of particle motion that is 
persistent relative to Brownian motion, and where a{t) 
decreases (Regime III); and a fourth regime where the 
non-Gaussian parameter has decayed back to very small 
values so that the particle motion is nearly Brownian 
(Regime IV). The four regimes are examinated in detail 
in the following subsections. 




FIG. 4. Classification of dynamic regimes. MFPT r(e) 
(solid curve) plotted versus e, and the non-Gaussian param- 
eter a{t) (dashed curve) plotted versus t, for T — 0.451. 
Inset: Comparison between r(e) (solid curve) and (r 2 (t)} 
(dot-dashed curve) on a linear scale. 

We emphasize that while some parallelism exists be- 
tween the mean square displacement and the first-passage 
time, the inset of Fig. ^ comparing these quantities shows 
that they are not equivalent. However, if the distribu- 
tion of particle displacments G s (r, t) were always exactly 
Gaussian (i.e. a(t) = for every t), then a simple in- 
verse function relation should hold between these quan- 
tities. G s (r, t) obeys a "scaling relation" if we can rescale 
G s (r,t) as, 



G 



(3.2) 



where / i s so me function. Simple dimensional analysis 
based on (|3.2|) implies, 



{r\t))^f 



(3.3) 



4 



and the scaling of the first-passage time with e, 

r(e) oc e 1/u . (3.4) 

In practice, these idealized scaling relations are re- 
stricted to certain time and space scales. Scaling with 
v = 1 is observed in the short time inertial regime. This 
result can be inferred from the general relation between 
the second moment of G s (r, t) and the particle velocity, 



(3.5) 



and the constancy of the velocity autocorrelation func- 
tion at short times [p|]43| , 



(3.6) 



(r(0) • r(i)) -> (r(0) • r(0)) = v 2 a t^O 



where v is the average particle velocity. Integrat- 
ing Eq. |3.5| over a short time interval obviously gives 
(r 2 (i)) ~ t 2 or "ballistic- like" motion. Actually, this scal- 
ing is just a consequence of the existence of equilibrium, 
and should not be construed as necessarily implying the 
absence of interparticle interactions at short timescales 
EH . A Gaussian form for the van Hove correlation func- 
tion in this short time regime is ensured by the Maxwell- 
Boltzmann distribution for the particle velocities. In the 
opposite extreme of very long times, the central limit the- 
orem governing the sum of independent particle displace- 
ments implies that the particle displacement distribution 
is Gaussian, and that v — 1/2. Transient scaling regimes 
can be observed at intermediate time scales, however. We 
refer to particle displacements as "persistent" relative to 
Brownian motion if v > 1/2, or "localized" relative to 
Brownian motion if v < 1/2. 



A. Inertial regime 

In the limit of very small e we probe the fast micro- 
scopic dynamics associated with the decorrelation of the 
particle momenta. It is difficult to probe this decorrela- 
tion directly using the first-passage time approximation 
to the dynamic entropy. Fig. |^ indicates that our ap- 
proximation for the dynamic entropy appears to diverge 
for e — > 0, so our approximation must break down in this 
limit. Gaspard and Wang have pointed out before that 
the first-passage time approximation to the dynamic en- 
tropy breaks down in this limit (l4| , so this shortcoming 
is to be expected. An estimate of the expected plateau in 
iS^e) corresponding to the Kolmogorov-Sinai entropy can 
be obtained by determining a cut-off time of the first- 
passage time distributions in the fast dynamics regime. 

The dynamics in this regime is examined by setting the 
magnitude of the first-passage sphere (e-sphere) about 
the center of each particle to be small enough that a 
collision does not usually occur before the particle leaves 
the e-sphere (see Fig. 1). By focusing on the particle first- 
passage time in this regime we identify a time and space 



scale over which particle velocities begin decorrelating. 
The first-passage time in this regime is insensitive to the 
type (A or B) of particle since they have the same mass. 
This is apparent in Fig. |^ where the first-passage time 
distributions for both the A and B particles are shown 
for the coldest run and e = 0.1. 
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FIG. 5. Probability of first-passage time Pt(t) for the A 
and B particles at T — 0.451 and e = 0.1. INSET: Same data 
represented in a semi-log scale. 
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FIG. 6. r versus e in the small e regime, 0.001 < e < 0.15, 
for T=0.468. Inset: log-log plot of the same data. Note 
that the asymptotic linear scaling of r(e) = 1.14e (solid lines) 
breaks down around a value e v » 0.05, corresponding to 
T v ~ 0.6. e v and t v are approximately independent of tem- 
perature (see also Fig. 3). 

The idealization of ballistic particle motion implies the 
time r(e) it takes for the particle to exit the sphere of ra- 
dius e equals e/vg. Fig |6| shows an expanded plot of r 
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versus e in the small e regime where this linear depen- 
dence of t on e is apparent. A non- linear dependence 
of r on e develops as the particle velocities decorrelate 
at a "velocity decorrelation scale" e v w 0.05. This dis- 
tance corresponds to the time t v w 0.6 (on the order of 
10~ 13 s in argon units). We find that e v and t v are ap- 
proximately independent of temperature in the present 
simulation (see also Fig. 3). 
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FIG. 7. First-passage time distributions Pe(t) in the iner- 
tial regime. Main figure contains six different temperatures 
at e = 0.1 plotted log-log. Solid line indicates Eq. 3.7 where 
the constant of proportionality has been adjusted to best fit 
the data. The long-dashed line indicates the cut-off distribu- 
tion function Eq. 3.8 with no free parameters. Inset: Pe(t) at 
T — 0.468 for e = 0.03. This value of e is less than the veloc- 



ity decorrelation scale e v indicated in Fig. 6. Eqs. 3.7 and 3.8 
are also shown in comparison with the simulation data. Note 
the absence of the long tail in the inset distribution. 

We obtain further insight into the decorrelation of par- 
ticle velocities in the inertial regime by examining the 
distribution of first-passage times, P e (t), as a function 
of e. Fig. 7 shows P e (t) for e = 0.1 and e = 0.03 for 
T = 0.468. We obtain a good approximation to P e (t) 
in the inertial regime through a Gaussian approximation 
for G s (r,t) |Q in conjunction with the first-passage dis- 
tribution for Brownian paths |^5| . The first-passage time 
distribution P e (t) scales as 



P e (t) ~ exp 



2(&)) 



(3.7) 



at short times and decays exponentially at long times, 
P e (t) ~ exp [— t/r(e)] J45|. We then introduce the ap- 
proximation 



A 



-(-)- — 

2\v tJ r(e) 



(3. 



where A is a normalization constant, (i? 2 ) 1 / 2 in Eq. 
has been replaced by v a t based on the assumption 



3.7 



that 



the particle displacements are nearly "ballistic" in the 
inertial regime. Note that the limiting expression for 
P e (t) given by Eq. 3.7 can also be deduced by assuming 
a Maxwell velocity distribution (one-dimensional due to 
the near rectilinear motio n) w ith the velocity replaced by 
e/t. Fig. 7 shows that Eq. (solid line) does not provide 
an accurate fit to the data using vo = 1/1.14 from the 
linear slope in Fig. ^, while the approximate expression 
Eq. 3.8 fits well for e less than the "velocity decorrela- 
tion scale" e v w 0.05. For e > e v , P e (t) shows evidence 
of developing a power-law "tail" at long times, as seen 
in the main part of Fig. 7. This tendency becomes more 
developed at larger e, as discussed in the next subsec- 
tion. The inset of Fig. 7 shows first-passage time data 
for e < e v where the tail is nearly absent and the cut-off 
is evident. At this point, we note that in liquids hxs 
characteristically scales as the inverse of the average in- 
terparticle collision time and thus has the interpretation 
of a microscopic "collision rate" JlJ,[l4| . The time t v can 
be considered as an average particle collision time so that 
an inverse relation betwen hxs and t v is expected. Wc 
do not consider the relation between t v and Uks, since 
simulations over a broader temperature range and the di- 
rect calculation of the Kolmogorov- Sinai entropy through 
Lyapunov exponent spectra are required for such 

an investigation. 



B. Particle localization regime 

The tendency of particle motion to become increas- 
ingly localized, as emphasized by the mode-coupling the- 
ory, is a conspicuous feature of experimental and sim- 
ulation data on supercooled liquids. The "plateau" in 
the mean square displacement log-log plots (see Fig. 2) 
indicates transient particle localization or "caging" , and 
the persistence of this plateau increases with decreasing 
T |4q . Next, we utilize the dynamic entropy concept to 
quantify this particle localization. 

An increase in the slope of log r(e) versus log e in Fig. || 
provides evidence for localization. A numerical differenti- 
ation of the data in Fig. || is shown in Fig. ||a where A(e) 
denotes the logarithmic derivative A(e) = dlogr/dloge. 
The scaling behavior in Eq. 3.4 implies that A(e) corre- 
sponds to the fractal dimension \jv of the particle tra- 
jectories. Fig. 8a shows that for all T, A(e) — » 1 for small 
e, and A(e) — ► l/i/ for large e. At intermediate values of 
e, we observe that A(e) develops a maximum at e c , cor- 
responding to the inflection point in Fig. |^. This length 
scale defines a distance that is difficult for the particle to 
exceed, and thus we define e c as the "cage" size. The inset 
in Fig. ||a shows that e c decreases with temperature, so 
that increased particle confinement occurs with cooling. 
Independent evidence indicating the significance of this 
characteristic scale is discussed later in this subsection. 

We denote the value of A(e) at e c by the "localiza- 
tion parameter" A(T) = maxA(e) = A(e c ), and its 



6 



T-dependence is shown in Fig ^a. The value of A(T) 
increases with cooling, consistent with increasing parti- 
cle localization {v < 1/2). The relatively noisy data in 
Fig. |b can be fitted by a power law, A(T) = 2 + 0.15(T- 
T c )-" 03 with T c = 0.435. 

Our identification of the cage size e c from the maxi- 
mum value of A(e) is further supported by the examina- 
tion of the first-passage time distribution P e (t) for e near 
e c . We find that P e (t) develops a long time power-law tail 
in this intermediate regime which is symptomatic of the 
development of intermittency in particle motion and par- 
ticle localization pq] . In Fig. ^ we show for T = 0.451, 
P e (t) at several values of e near e c = 0.21 ± 0.02. The 
apparent power law for the P e (t) tail varies with e and 
has a value near 2 for e « e r . 
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FIG. 8. (a) A(e) = dlogr(e)/dloge) for (from top to bot- 
tom) T = 0.451,0.457,0.468,0.480,0.505, and 0.550. Inset: 
"Cage" size e c versus T. (b) A = max A(e) plotted versus T. 
The dashed line is a power law fit to the data as indicated in 
the text. 

This power-law tail behavior in P £c (i) is shared by the 
four coldest temperatures, as shown in Fig. 10, where 
the asymptotic behavior of the distributions is seen to 
be numerically very similar. The difference in their first 
moments [i.e. r(e c ) in Fig. 3] reflects both the differences 
between the distributions in Fig. 10 at short times, and 
the asymptotic cut-off in the distributions that is difficult 
to resolve numerically. 
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FIG. 10. First-passage time distributions at the cage scale, 
e — e c , for four different T. 

Odagaki |3(| has suggested that the second and first 
moments of the first-passage time distribution at the 
scale of one interparticle distance diverge at T c and at 
the glass transition temperature T gi respectively, and a 
recent paper by Hiwatari and Muranaka fill supports 
this glass transition scenario. A transition to intermit- 
tent particle motion at the glass transition has also been 
suggested by Douglas and Hubbard |37f| . Our observa- 
tions are consistent with the growth of intermittency of 
particle motion as T c is approached, but we cannot con- 
firm theoretical predictions of a dynamical transition in 
the degree of intermittancy until lower temperatures are 
examined. An accurate test of these predictions will re- 
quire carefully equilibrated data below T c , beyond the 
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temperature range of the simulations analyzed here. We 
point out that, in the present system, the scale at which 
this intermittency occurs is substantially smaller than 
one interparticle spacing, and instead pertains to motion 
at the scale of the cage size, e c . 

C. Regime of persistent particle motion 

Previous work has shown that particle motion becomes 
increasingly collective in this supercooled liquid and that 
an important mode of motion at the scale of the interpar- 
ticle distance involves the string-like collective motion of 
particles. A visualization of this process |^0| suggested 
to us that this process becomes increasingly "coherent" 
or "jump-like" at lower temperatures and this tendency 
towards "coherent jumping" has been noticed in a num- 
ber of other physical systems (melting of hard disks , 
hexatic liquids ||52| and ordering plasmas (53)). In this 
subsection we utilize S(e) to further quantify this effect. 

Regimes I and II were defined by characteristic spatial 
scales at which changes occur in the first-passage time 
distributions. However, the long run times of the simula- 
tions analyzed here necessitated the storing of configura- 
tions on a logarithmic, rather than linear, time scale, for 
all but the coldest simulation || . As a consequence, first- 
passage time distributions cannot be obtained over a con- 
tinuous range of e for large e. In the absence of this infor- 
mation, we roughly identify Regime III by the tendency 
for the non-Gaussian parameter a(t) to decrease (see 
Fig. 4). In Fig. 11 we show the apparent fractal dimen- 
sion A(e) in Regime III for the highest and lowest temper- 
atures. (The results for all T are shown over an extended 
scale in Fig. 8). Notice that persistent particle motion 
(A(e) < 2) develops for e > 0.6 in Fig. 8 (although A 
remains near 2 at the highest temperatures) and conse- 
quently we show A(e) for the range 0.7 < e < 1.0 in 
Fig. 11. The value of e at which A(e) = 2 provides a 
more precise estimate for the beginning of Regime III. 
Although the data is noisy, we find that A(e) w 2 within 
numerical error, and is nearly independent of e for high 
T. Particle displacement at high temperature is then 
reasonably approximated by Brownian motion on these 
spatial scales. However, a substantially smaller average 
value of A(e) 1.7 (dashed line in Fig. 11) is found for 
the lowest T. Thus, we find the particle motion becomes 
increasingly persistent on cooling. 

Why is persistent motion not observed in the mean 
square displacement? The tendency for persistent par- 
ticle motion is not apparent in (r 2 (t)) shown in Fig. ^ 
or in the inset of Fig. ^| because, when averaging over 
the squared displacements, the contribution of the few 
particles that are at any given time moving persistently 
is "washed out". However, these particles give a large 
contribution to the mean first-passage time, so that this 
quantity is therefore a sensitive indictor of persistent par- 
ticle motion. 
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FIG. 11. Apparent fractal dimension A(e) of particle dis- 
placements in Regime III. From Fig. 8, highest and lowest 
simulation temperature T = 0.5510 and T — 0.4510 are shown 
for 0.7 < e < 1.0. 

We can obtain insight into the emergence of persis- 
tent particle motion in our cooled liquid from idealized 
models of Brownian motion subject to potential fluctu- 
ations. If the potential fluctuations are quenched, there 
is a tendency towards particle localization Q, but the 
occurrence of fluctuations in the potential in both space 
and time can lead to persistent particle motion. In par- 
ticular, if the fluctuations are delta-correlated in both 
space and time, then the exponent A equals 3/2 [ p5[ . 
The effects observed in Ref. (HHj are qualitatively consis- 
tent with our understanding of the origin of correlated 
motion. At short times, the existence of relatively im- 
mobile particles leads to a randomly fluctuating field felt 
by those particles free to move at a given point in time. 
This spatially fluctuating field is responsible for the par- 
ticle caging or localization on timescales short compared 
to the decorrelation time r Q of the "structural fluctu- 
ations" (associated with the relatively immobile parti- 
cles). At longer times, the formerly immobile particles 
become mobile and the potential field fluctuates in time, 
leading to an enhancement in the particle displacement. 
At still longer times, thermal fluctuations restore equilib- 
rium and particle displacement ultimately becomes dif- 
fusive. 

Finally, we point out that particle motion can be per- 
sistent even in the absence of a secondary (so-called "hop- 
ping") peak in G s (r,t). Instead, persistent particle mo- 
tion contributes to a long tail in G s (r, t) in the temper- 
ature range of the present system q^|> This long tail 
sharpens up and becomes a secondary peak at lower T 
p6| , indicating the increased contribution of collective 
particle motion to transport below T c . 
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D. Large scale particle displacement 

Particle displacement in a liquid at large scales is de- 
scribed by Brownian motion so that 5(e) should scale 
asymptotically as e 2 for large e, regardless of temper- 
ature. It is apparent in Fig. g that the data in the 
asymptotic diffusive regime is limited, especially at lower 
temperatures. Previous work has shown that there is a 
tendency for "mobile" particles, which dominate trans- 
port in cooled liquids near T c , to move an interparticle 
distance during the time in which they are mobile 
This happens because the particles tend to move between 
local minima in the potential surface describing the in- 
terparticle interaction J57],Q. This feature is especially 
apparent in the string-like particle motion noted before, 
where it has been observed that the strings "disintegrate" 
once an interparticle displacement has been achieved Q . 
Thus, one interparticle distance is taken to be the mini- 
mal scale of the large scale particle displacement regime 
(Regime IV). 
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FIG. 12. Dynamic entropy at the scale of one interparticle 
separation, S(e = 1), versus T — T c , with T c = 0.435. Di- 
amonds refer to the simulation data for S(e = 1), and the 
dashed lines refer to the power law, S(T) ~ (T — T c ) 2 5 . We 
also observe this scaling for S(e — 0.4) (circles). The in- 
verse of the structural relaxation l/r a (triangles) follows a 
power law l/r a ~ (T — T c ) 2 ' 8 . Inset: Schematic indication of 
expected temperature variation of dynamic entropy Hks at 
small scales. Note that the extrapolated ergodic-nonergodic 
transition in the e — > limit should occur at a temperature 
T < T c . 

In the previous discussion, we have established that 
5(e) is insensitive to temperature for small e over the tem- 
perature range investigated. This insensitivity accords 
with the expected variation of dynamic entropy hxs- 
At higher T we expect the dynamic entropy to saturate, 
while a decrease should accompany the more restricted 
particle motion at lower T [p5[ . A variation similar to 
that shown in the inset of Fig. 12 has been established 
for the ordering of the XY model [Q . The investigation 
of the anticipated ergodic-nonergodic transition and its 
possible relation to a vanishing of hxs requires the calcu- 
lation of the full Lyapunov spectrum, which is currently 
prohibitive for a system of the present size. However, we 
can investigate 5(e) at a larger scale on the order of an 
interparticle spacing. This should be interesting because 
of the strong T-dependence in 5(e) at this length scale 
noted above (see Fig. 3). 

In Fig. [l^ we examine the T-dependence of 5(T; e) 
at the scale of one interparticle separation e = 1. Wc 
observe that 5(T;e = 1) obeys a power law to a good 
approximation over the temperature range investigated, 
and that 5(T; e = 1) seems to extrapolate to zero at the 
mode-coupling temperature T c = 0.435. As shown in 
Fig. 12, a reasonable fit to the data is obtained with the 
relation 



5(e = 1) ~ (T - T c 



\2.5 



T c = 0.435. 



(3.9) 



The scaling of 5(r; e = 1) is compared in Fig. [T^ to the 
structural relaxation time r a describing the decay of the 
intermediate scattering function. Although the scaling 
of the two quantities is qualitatively similar, the best fit 
exponent for T a has the somewhat larger value of —2.8. 

The power law scaling of 5(e = 1) with temperature 
is not obvious since if the particle displacement were ex- 
actly described by Brownian motion, then r(e) would 
scale in inverse proportion to the diffusion coefficient. A 
determination of the diffusion coefficient for A particles 
is obtained by a simple least-squares fit (not shown) of 
the long time data to the function (r 2 (t)) = A + 6Dt. 
This fitting gives 



D~{T- T c f 



(3.10) 



A more refined estimate by Kob and Anderson [ p9[ on 
a smaller (1000 particles) system gave an exponent 2.0 
for the A particles. The diffusion data thus scales with a 
fractional power of the structural relaxation time, 
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D 



-(2.1/2.8) 



-0.75 



(3.11) 



over the temperature range investigated. Since evidence 
supports a common temperature scaling of r a and the 
fluid viscosity i] f39f , the observation implied by Eq. 
(3.11) is consistent with the breakdown of the Stokes- 
Einstein relation in real and simulated supercooled liq- 
uids pp|j6l[ . We therefore conclude that 5(e — 1) scales 
neither exactly like the inverse structural relaxation time 
1/t q nor the diffusion coefficient D. Other characteristic 
times exist for this liquid. It has recently been reported 
for the same simulations investigated here that the time 
scale on which particle displacements are most correlated 
scales as a power law with (T — 0.435), with an exponent 
7 = 2.3 ± 0.2 0. This exponent notabl y a grees within 
numerical error with the exponent in Eq. We further 
note that the time t* at which a(t) is a maximum (see 
Fig. 4) appears to diverge at T c with an exponent 1.7 jjj. 

Since the temperature dependence of 5(e) is largest in 
Fig. H for e — 1 and very small for e — > 0, it is natu- 
ral to consider the temperature dependence of 5(e) for 
many fixed e to determine how this crossover occurs. In 
Fig. [l| we see that the scaling 5(T) oc (T — T c ) 2 5 holds 
for e greater than the cage size e c and sufficiently small 
reduced temperature. We also find that for any temper- 
ature, the scaling fails for e < e c . This provides another 
method for determining the cage size. 
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IV. CONCLUSION 

We have studied the length-scale dependence of the dy- 
namic entropy 5(e) in a molecular dynamics simulation 
of a model supercooled binary Lennard- Jones liquid. The 
simulations were performed above both the glass transi- 



tion temperature and the mode-coupling critical temper- 
ature and correspond to equilibrated liquid states [|3j . 

5(e) as estimated by the MFPT provides a tool for 
identifying characteristic length and time scales of the 
dynamics of liquids and provides a means for quantifying 
the degree of correlated motion occurring in supercooled 
liquids. At very small e we observe a decorrelation of 
particle velocity and an 5(e) which is insensitive to tem- 
perature and e variation. A decorrelation time associated 
with an average interparticle collision time is identified. 
At intermediate values of e we observe a sharp drop in 
5(e) with increasing e, indicating that the path motion 
is more "stochastic" at these length scales. 5(e) is found 
to have a strong temperature dependence as well in this 
regime. The logarithmic derivative A = d log r/d log e 
becomes a maximum at a characteristic e value that is 
identified with the particle "cage" size e c , since parti- 
cle localization is maximal at this point. The scaling of 
5(e) with temperature at fixed e gives an independent 
confirmation of our estimation of e c because it exhibits 
a qualitatively different dependence on temperature for 
e > e c and e < e c . The localization parameter A = A(e c ) 
increases and the cage size e c decreases as the liquid is 
cooled. The distribution functions for the first-passage 
time at the scale of the cage size, P £c (t), in the coldest 
runs exhibit a long power-law tail, consistent with sugges- 
tions that there is growing intermittency in the particle 
displacements in glass- forming liquids ^,^]. This fea- 
ture requires further study in cooler liquids to check the 
predictions of these models. At still larger scales (still 
less than one interparticle separation) we observe persis- 
tent motion which follows the transient particle "caging" . 
5(e) obeys a power law 5(e) ~ e~ 1/V for e in the range 
(0.7, 1) where the apparent fractal dimension of the par- 
ticle trajectories A(e) = 1/v ranges from w 2 to w 1.7 
as the temperature is lowered. Thus, we observe a ten- 
dency for the particle motion to acquire a persistent char- 
acter relative to Brownian motion as it is cooled. This 
is consistent with the previous observation of correlated 
string- like motion in this liquid 

General arguments suggest that the dynamic entropy 
at microscopic scales decreases at low temperatures, but 
a slower variation should be obtained at higher temper- 
atures, as in the present molecular dynamics calculation 
(see Fig. 2). The situation is not so clear for the first- 
passage time dynamic entropy 5(e) at the scale of one 
interparticle separation e — 1. In this case 5(e = 1) is 
inverse to the "average interparticle exchange time [p2[ " , 
r(e = 1). We find that 5(e) vanishes as a power law, 
5(e = 1) oc (T — T c ) 2 - 5 when the mode-coupling tem- 
perature T c is approached from above. Thus, the cooled 
liquid has the appearance of approaching an ergodic to 
non-ergodic transition as T —> T c , when viewed at the 
scale e = 1. This is consistent with the predictions of 
mode-coupling theory. 

The first-passage time can also be utilized to obtain in- 
formation about the spatial dependence of mobility fluc- 
tuations in cooled liquids. Perera and Harrowell |63[, for 
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example, have examined the position dependence of r at 
the scale of one interparticle spacing in a two-dimensional 
soft-sphere supercooled liquid and found a tendency for 
particles of relatively high and low "mobility" (i.e., small 
and large r, respectively) to cluster as T is lowered. A de- 
tailed study of spatial correlations of first-passage times 
in the present system will be presented elsewhere J64j . 

Future work should examine our approximate expres- 
sion for S(e) in the e — > limit through independent 
calculation of the Kolmogorov-Sinai dynamic entropy, 
S(e — > 0) = Iiks- The temperature dependence of h^s 
in cooler liquids (T < T c ) should be examined to de- 
termine if there is a tendency for the "bare" dynamic 
entropy to vanish at a lower glass transition tempera- 
ture (see inset of Fig. 12). Recent work has established a 
phenomenological relation between dynamic entropy has 
and the equilibrium entropy in simulations of model liq- 
uids at relatively elevated temperatures Ji6} |. A decrease 
in hj^g at lower T should be accompanied by the devel- 
opment of collective motion at short times. Such motion 
has been reported in Ref. jl9|. We expect this change 
in the short time dynamics to be relevant for interpret- 
ing the Boson peak phenomenon in cooled liquids J65| . 
Simulations have already shown that the fraction of un- 
stable modes f u in a cooled liquid decreases in parallel 
with hxs P5|] , and the vanishing of f u has been identi- 
fied with the temperature where the diffusion coefficient 
D vanishes f36|,|67|. 

Finally, we note that the calculation for S(e) can be 
extended to other dynamical variables associated with 
other transport properties (viscosity, thermal conductiv- 
ity, etc.) and these calculations should provide useful 
estimates of other characteristic space and time scale in 
cooled liquids |Q . We emphasize that although the def- 
inition of S (e) is motivated by dynamical systems theory 
concepts, this quantity defines an independently interest- 
ing measure of correlated motion in liquids that does not 
rely on the approximation relating S(e) to h(e). 

Corresponding author: sharon.glotzer@nist.gov. 
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